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ABSTRACT 



> 

We explore the intergalactic medium (IGM) as a potential source of the un- 
resolved soft X-ray background (XRB) and the feasibility to extract the IGM 
state and evolution from XRB observations. We build two analytical models, 
the continuum field model and the halo model, to calculate the IGM XRB mean 
flux, angular auto correlation and cross correlation with galaxies. Our results 
suggest that the IGM may contribute a significant fraction to the unresolved 
soft XRB flux and correlations. We calibrated non-Gaussian errors estimated 



against our 512 3 moving mesh hydro simulation and estimate that the ROSAT all 
sky survey plus Sloan galaxy photometric redshift survey would allow a ~ 10% 



accuracy in the IGM XRB-galaxy cross correlation power spectrum measure- 
ment for 800 < / < 5000 and a ~ 20% accuracy in the redshift resolved X-ray 
emissivity-galaxy cross correlation power spectrum measurement for z < 0.5. At 
small scales, non-gravitational heating, e.g. feedback, dominates over gravity and 
leaves unique signatures in the IGM XRB, which allows a comparable accuracy 
in the measurement of the amount of non-gravitational heating and the length 
scales where non-gravitational energy balances gravity. 

Subject headings: Cosmology-theory: X-ray background, intergalactic medium, 
large scale structure 
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1. Introduction 

A large fraction of the intergalactic medium (IGM) is hot. The cosmic virial theorem 
predicts that the mean IGM temperature > 0.2 keV (Pen 1999). It emits X-ray through 
thermal bremsstrahlung and contributes to the soft (0.5-2 keV) X-ray background (XRB). 
About 80-90% of the soft XRB has been resolved into point objects such as AGNs (Hasinger 
et al. 1993). Various sources such as nearby low luminosity AGNs (Halderson et al. 2001), 
unresolved galactic stars (Kuntz & Snowden 2001), galactic gas, X-ray sources in external 
galaxies and the IGM (Pen 1999; Wu k Xue 2001; Croft et al. 2001) may contribute a 
significant fraction to the remaining 10-20%. To distinguish those possible components, one 
can combine the XRB mean flux, auto correlation (Soltan & Hasinger 1994; Sliwa, Soltan & 
Freyberg 2001) and cross correlation with galaxies (Almaini et al. 1997; Refregier, Helfand 
& Mcmahon 1997). In §2, we will build analytical models to estimate the IGM contribution 
to the XRB flux and correlations. 

In pure gravitational clustering where the only source of thermal energy is shock heat- 
ing from collapse, simulations show that the gas correlation function £ gas follows the dark 
matter correlation function down to scales where the gas over-density correlation £ gas ~ 10 3 
(Pen 1999). Non-gravitational heating (feedback) can rearrange the gas distribution on small 
scales. By searching for the scale where feedback dominates over gravity, one can robustly 
measure the strength and history of non-gravitational heating. The IGM XRB is sensitive 
to small scale gas structures due to the X-ray emissivity dependence on density squared. 
This makes the IGM XRB a potentially powerful probe for studying feedback from galaxies 
on the IGM. From the estimation of the IGM XRB flux, Pen (1999); Wu, Fabian & Nulsen 
(2001); Voit & Bryan (2001) suggested that, if only gravitational heating exists, the soft 
X-ray emission that the IGM produces would exceed the observational limit. To suppress 
the IGM clumping and reduce the X-ray emission, a significant amount of non-gravitational 
injection energy ~ 1 keV/nucleon is required. Croft et al. (2001); Dave et al. (2001); Phillips 
et al. (2001) argued from simulations that such energy injection is not necessary. In their 
simulations, a large fraction of baryons (about 70%) are too cold to contribute significantly to 
the XRB. Despite the stability problem of cool gas, these findings may indicate a degeneracy 
between non-gravitational heating and the fraction of IGM contributing to the XRB. The 
IGM XRB auto correlation function (ACF) and cross correlation function (CCF) with galax- 
ies have different dependences on the IGM thermal state and are capable of breaking this 
degeneracy. Furthermore, with galaxy photometric redshift data, the redshift resolved IGM 
X-ray emissivity-galaxy cross correlation and emissivity auto correlation can be extracted. 
This tells the evolution of the IGM state. In §3, we will forecast the sensitivity of ROSAT 
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all sky survey (Voges et al. 1999) and SDSS 3 and test the feasibility to constrain the IGM 
thermal history from correlations. This estimation depends on the XRB non-Gaussianity 
since it affects the error analysis. We analyze our 512 3 moving mesh hydro simulation to 
quantify this effect. 



2. Analytical models for the IGM XRB 

The X-ray emissivity in the band E x < his < E 2 is given by: I x ^ 2AxlO- 27 c z T 1 / 2 [ex.p(-E 1 /T)- 
exp(-£ 2 /T)]n 2 = Ag(T)n 2 erg cm~ 3 s' 1 (Tucker 1975). Ag(T) is the cooling coefficient 
in the energy band E 1 -E 2 . A mean gaunt factor 1.2 is adopted. c z ~ Z/Z Q (4keV /T) + 1 
is a fitting formula to result of metal cooling from Raymond, Cox & Smith (1976). We will 
use Z = \Zq as the lower plausible limit on the IGM inferred from clusters of galaxies. We 
choose a flat ACDM cosmology with Q = 0.37, Q A = 0.63, Q B = 0.04, a 8 = 0.9, h = 0.7 
and n — 1 to calculate the soft XRB statistics. Hereafter, we always consider the comoving 
emissivity 1% — I x /(1 + z ) 3 - 

The IGM XRB flux F x (0) = J /£(x#)/[4tt(1 + zf]dx- X is the comoving distance. The 
mean XRB flux F x = J /£/[4tt(1 + z) 2 ]d X - The galaxy surface density £(0) = / g[l + 
$G(x9)]dz. is the galaxy redshift distribution and 5q is the galaxy number over-density. 
The mean surface density £ = f ^dz. We define the dimensionless fluctuations = 
F/F - 1 and A E = £/£ -I. The IGM XRB ACF and CCF with galaxies are defined by 
w x {0) = (A F (9 1 )A F (§ 2 )) and w XG (9) = (A F (9 1 )A^(9 2 )), respectively. Here, 9 1 -9 2 = cos6». 
C x (l) and C XG (l) are the corresponding angular power spectra, respectively. 

These 2D correlations are determined by the corresponding 3D correlations such as the 
emissivity ACF £x( r ) = (^x(x)5x(x + r)) and emissivity-galaxy CCF ^ G (r) = (5 x (x.)5g(x + r)) 
or their corresponding power spectra P x (k,z) and P X c{k,z). 5 X = I x /I x — 1. At small 
angular scales we use the Limber's equation to obtain 

Here, Xre is the comoving distance to the reionization epoch. [z±, z 2 \ is the redshift range of 
the galaxy survey adopted. 



3 SDSS, http://www.sdss.org/ 
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For these statistics, we can treat the IGM either as a continuum field with a density 
and temperature distribution or as discrete halos. From these two viewpoints, we build two 
analytical models: the continuum field model (§2.1) and the halo model (§2.2). 



2.1. The continuum field model 

For the [0.5, 2] keV X-ray band, we can approximate the gas temperature T ~ 1 keV 
due to the following arguments. (1) when T <C 1 keV, A^(T) decreases exponentially and 
we would expect that too cold gas does not contribute significantly to the soft XRB. (2) 
When T ^> 1 keV, A^(T) drops as T" 1 / 2 . Furthermore, gas with T ^> 1 keV is rare. So 
the contribution from gas with T ^> 1 keV to the soft XRB is small. (3) A^?(T) peaks 
at T ~ 1 keV. The density weighted temperature is ~ 0.4 keV (Zhang & Pen 2001) and 
the density square weighted (roughly emissivity weighted) temperature ~ 1 keV. So, we 
expect most contributions to the soft XRB to be from gas with T ~ 1 keV. Around this 
temperature, the emissivity has only a weak temperature dependence. So, we approximate 
T = 1 keV and define a X-ray weighted gas fraction f^ s = ^Ix/[A^?(1 keV)n 2 ]. From 
the above argument we expect f* as ~ 1, but throughout this letter, we treat as a free 
parameter to be determined by observations. At large scales, gas generally follows dark 
matter. At smaller scales, the gas density fluctuation is suppressed by the gas pressure. 
This effect can be modeled through a window function W g (r) such that the gas over-density 
is a convolution of the dark matter density and W g (r). The FWHM of W g (r) is the scale 
where non-gravitational processes such as feedback begins to dominate over gravity. In 
Fourier space, the Fourier component of the gas over-density S gas (k) = 5{k)W gas {k) where 
5(k) is the Fourier component of the dark matter over-density. A Gaussian window function 
W gas (k) = exp(— fc 2 //Cg as ) is choose from the literature (Gnedin & Hui 1998). With the 
correlation between gas and dark matter, one can calculate the XRB statistics using the 
knowledge of dark matter density field. 

We apply the extension of the hierarchical model (Fry 1984) in the highly non-linear 
regime: the hyper-extended perturbation theory (HEPT) (Scoccimarro & Frieman 1999) to 
calculate P x and Pxg- 

P x (k, z) ~ 4 1 / B 4 (k u k 2 , k 3 , k 4 ; z)W 1 W 2 W 3 W 4 d 3 k 2 d 3 h (3) 

Pxc(k, z) ~ - 1 / 65 3 (ki, k 2 , k; z)W 1 W 2 d 3 k 1 . (4) 

Here, Wi = W ga , s (ki), ki+k 2 = — k 3 — k 4 = k. The gas clumping factor o"g as = J P gas (k)k 2 dk/ (2n 2 ) 
with Pgas(^) as the gas density power spectrum. The bispectrum B 3 oc P 2 (k) and polyspec- 




Fig. 1. — The IGM cumulative contribution to the XRB. Contribution with different k t 
looks similar. 
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Fig. 2. — IGM XRB angular ACF w x (9) and CCF w X g(Q) calculated from the contin- 
uum field model (left panel) and the halo model (right panel). The solid, dot, dash lines 
correspond to k gas = 8, 4, 3/i/Mpc ( corresponds to non-gravitational heating energy per nu- 
cleon £ NG ~ 0.3,0.6,0.9 keV) in the left panel and r c {M 8 ,z = 0) = 0.5, 0.75, 1.0/r 1 Mpc 
(E NG ~ 0.5, 0.6, 0.8 keV) in the right panel. k gas and r c are the scales below which feedback 
dominates over gravity. The slopes of w(8) predicted by the halo model are steeper than the 
continuum field results. This difference is likely due to the weak correlation between density 
and temperature, which is omitted in the continuum field model. 
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trum, the dominant term in the expression of Px(k, z), oc P 3 (k) terms are calculated from 
HEPT. Here, P(k) is the dark matter density power spectrum. See Zhang & Pen (2001) for 
more detailed explanation. We consider a flux limited galaxy survey SDSS, take the SDSS 
galaxy distribution as ^ = 2 ( z /i 412) 3 ex P[~ (1-412^/zm) 3 / 2 ] (Baugh & Efstathiou 1993) and 
adopt z m = 0.43 as the median redshift of Sloan galaxy photometric redshift distribution 
(Dodelson et al. 2001). We assume that k gas does not change with redshift and a constant 
bias model 5q = b5 for galaxies. 



2.2. The halo model 

The gas profile in a halo and halo mass-temperature relation determines the X-ray lumi- 
nosity. The halo mass function and halo-halo correlation enable us to calculate their collective 
effects to the XRB. Similar methods have been applied to the dark matter correlation (Ma 
& Fry 2000) and the Sunyaev-Zel'dovich effect (Komatsu & Kitayama 1999). 

-1 
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The gas core ra- 



We adopt the electron number density profile n e = no 
dius r c is analogous to k gas and corresponds to the scale below which feedback significantly 
changes gravitational clustering. We assume that gas accounts for f^/f2 = 11% of the 

halo mass. The gas temperature is determined by the virial theorem through the relation 

3/2 



M 
M 8 



r/kcv 



, ., o „ v , (Pen 1998). M 8 is the mean mass contained in a 8h 1 Mpc 

sphere today. The distribution of halo comoving number density n as a function of halo 
mass M and z is given by (Press & Schechter 1974): = (f ) 1/2 {^^\^\ exp(-Jj). 
Here po is the present mean matter density of the universe. <r(M, z) is the linear theory 
rms density fluctuation in a sphere containing mass M at redshift z. 5 C is the linearly ex- 
trapolated over-density at which an object virializes. Its dependence on cosmology is weak 
and we will adopt the value 1.686, which is the 5 C for a Qq — 1-0 CDM universe Eke, Cole 
& Frenk (1996). Mo & White (1996) related the halo-halo correlation P C (M U M 2 ) with the 
underlying dark matter correlation P(k) by a linear bias: P c (k, M 1? M 2 ) = P(k)b(M 1 )b(M 2 ). 
We adopt the NFW profile (Navarro et al. 1996) with a compact parameter c = 5 to calculate 
the XRB-dark matter cross correlation. 

In this model, the variance of the gas density a 2 gas oc n^r 3 [tan -1 (r A / r c ) — p ] j^dM. 
ta is the virial radius and we calculate it using the fitting formula of Eke, Cole & Frenk 
(1996). Defining <5 x (k, M) and <5(k, M) as the Fourier transform of I%(r, M) /I$ and 5(r, M) = 
Pdm{ y -i M)/ Pdm of each halo, respectively, we obtain 
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nil 

%(k) m dM + P(k) 



Jw low 



(5) 



Pxs(k) = 



p{k) 



(6) 



ow 
oo 



nil 

S x (k)^dMb(M) 



5{k)^dMb{M) 



The integrals in these equations depend strongly on the halo lower mass limit M\ ow , which 
can not be arbitrarily chosen. A smaller Mi ow will produce a bigger gas clumping factor 
cTg as since more gas contributes. It will also produce a smaller Px since gas in less massive 
halos is more diffuse and contributes a smaller fraction to the correlation than to the mean 
flux. This behavior contradicts the P x oc <r| as behavior we should expect. Thus M low must 
be determined independently. Here we adopt the model of Pen (1999). In this model, the 
energy injection from feedback such as supernovae winds expand the gas and produces a core 
(radius ri(M)) with a constant entropy. Then r^M) = r A (M) sets the value of M low . In this 
model, r 1 ~ 0.5r c oc T _1 oc M~ 2 / 3 . We further assume r c (M 8 ,z) = r c (M 8 ,z = 0)/(l + z), 
which corresponds to a redshift independent /c gas . 



2.3. Predictions 



Two models give consistent predictions for the XRB. More than 80% of the contri- 
bution to the soft XRB flux is from the IGM at z < 1 (fig. 1). The mean X-ray flux 

2 

F x c± (^)ICT 12 erg s~ x cm~ 2 deg~ 2 , which accounts for a significant, if not dominant, 
fraction of the unresolved XRB. The IGM XRB is homogeneous with a large amplitude 
(fig. 2) and is sufficient to explain the observed XRB correlations. The point-to-point 
XRB variance W x (0) ~ 10 and W XG (0) ~ 0.3. For 9 > l', we find that W XG oc 6 a 
(a ~ —1.1). The shape and amplitude of these properties have different dependences on 
gas parameters. (1) For the shapes, the larger A; gas (corresponds to a smaller r c and a 
larger a 2 as ), the steeper the correlations. (2) For the amplitudes, F x oc (/g^ s ) 2 o"| as while 
Px(k), Pxc(k), wx(0), wxg(Q) oc (/g^ s )°o"| as - Combining correlation data and mean flux 
data, one can distinguish from <r| as . With redshift resolved Px(k, z) or P X a(k, z), as can 
be obtained from XRB observations and galaxy surveys, one can infer cr 2 as (z) and then k gas (z) 
and r c (z) from their relations with <r| as (z). The interpretation of these quantities depends on 
processes changing the gas-dark matter correlation such as non-gravitational heating and pos- 
sibly gravitational shock-heating. Current simulations have difficulties to resolve these small 
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scales where these processes become dominant due to limited simulation mass and spatial 
resolution. Simulation results about the IGM XRB flux, hot gas fraction, etc. have not con- 
verged, so hereafter we focus on the estimation of the effect of feedback through our analytical 
models and use our hydro simulation to calibrate some XRB statistics. Feedback increases 
the gas temperature and thus /^ s , the fraction of gas contributing to XRB. But it dilutes 
the gas and reduces <7g as , resulting in a larger r c or a smaller k gas . From Pen's model, for a 
ACDM universe, we can estimate E^q(z), the non-gravitational energy injection per nucleon 
in unit of keV, from the relation <r| as ~ 90£l [2y/E NG exp(— E^g) + V^^c(y/E^ G )]/E^G). 
k gas (z) (or r c ) and E N q(z) then tell us the scales at which feedback dominates over gravity 
and the feedback strength. 



3. Extracting the IGM state and evolution 

X-ray sources and differential extinction of our galaxy make the measurement of the 
IGM XRB flux and ACF difficult. But the measurement of IGM XRB-galaxy CCF is much 
more robust. The direct observable in the CCF measurement is (F x (9i)T, G (9 2 )). We need 
to estimate the IGM contribution to the above property in order to infer w I x G A (9). (1) Our 
calculation suggests that the IGM XRB is sufficient to explain the unresolved X-ray flux 
and the cross correlation with galaxies. (2) Distant AGNs and galactic X-ray sources have 
almost no correlation with nearby galaxies with z < 1. (3) The CCF between galaxies 
and X-ray sources in extragalactic galaxies or nearby low luminosity AGNs is of the same 
amplitude of galaxy surface density ACF, which is one order of magnitude lower than the 
IGM XRB-galaxy CCF. So, even if they contribute a comparable amount to the XRB flux 
as the IGM, their contribution to the cross correlation is negligible. (4) For a low matter 
density universe, the CCF caused by the weak lensing of low redshift large scale structures to 
unresolved high redshift AGNs accounts for at most 1/3 of the correlation (Cooray 1999). In 
principle, combining the XRB mean flux, auto correlation and cross correlation measurement, 
the IGM XRB cross correlation can be determined. We neglect possible systematic errors in 
such measurements and take the ROSAT all sky survey and SDSS as our targets to estimate 
the statistical error in the IGM XRB CCF measurement. 

The ROSAT all sky survey (RASS) covers the whole sky (f^ ss ~ 1) in the 0.1-2.4 soft 
X-ray band in t — 1.03 x 10 7 s. At the energy bands below 0.5 keV, the galactic emission 
and absorption are strong and at bands > 1 keV, only a fraction of the sky can be used. So, 
we choose the band 0.5-1 keV for our analysis. A factor of 4 decrease in the photon count 
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Fig. 3.— A XRB map in our 512 3 MMH simulation(Zhang, Pen & Wang 2002). The 
parameters we adopted in this simulation are flo = 0.37, Q\ = 0.63, Qb = 0.05, h = 0.7, 
erg = 1-0, the primordial power spectrum index n = 1, box size L = 100/i _1 Mpc and smallest 
grid spacing 40/i _1 kpc. The XRB map integrates contributions up to z = 3 and has a width 
of 1.32 degree. The XRB flux is in units of erg s _1 cm -2 sr _1 . 




Fig. 4. — The XRB-galaxy cross correlation non-Gaussian corrections r)(l) from sample 
variance errors measured from our 512 3 MMH simulation. For Gaussian signals, rj(l) = 0. 
We make 40 XRB maps and 40 maps of the intergalactic gas surface density at the same 
fields of view to measure the XRB-galaxy cross correlation assuming that galaxy number 
density traces the gas distribution. We then measure the fluctuation of the cross correlation 
power spectra of these maps to obtain the XRB-galaxy cross correlation power spectrum 
non-Gaussianity r)(l). We find a strong non-Gaussian correction at scales 500 < I < 10 4 . 
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Fig. 5. — Statistical errors in the forecast ROSAT+SDSS measurement of - ^ C XG ( I) _Fjp M 
(left panel) and extracted Pxc{k,z) (right panel). In the left panel, error bars are for the 
central line with r c (Mg, z = 0) = 0.75/i _1 Mpc and we normalize its _Fjp M = 1. Solid 
error bars are Gaussian errors and the dot error bars are non-Gaussian errors. The effect of 
non-Gaussianity is dominant at large angular scales and is negligible at small angular scales 
where the dominant error is from noise. The redshift averaged gas parameter r c can be 
constrained with ~ 10% accuracy at la. The amount of feedback and the scale below which 
feedback dominates over gravity can be modeled with a comparable accuracy as Pxc{k, z). 
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rate is expected due to the choice of this band instead of the whole RASS band 4 . One full 
field view of RASS has a ~ 1 degree radius and only within the central 10-15 arc-min radius 
is the resolution better than 3 arc-min. Since we need high resolution to probe the IGM 
state, we will focus on these central regions. Since the separation of two successive scans 
is smaller than 10 arc-min, these high resolution regions cover almost the whole sky. But 
the number of photons received in these regions is about l/36-l/16th of the total number 
of photons received. We adopt the conservative fraction 1/36 for our following estimation. 
The SDSS covers / s |° ss ~ 1/4 fraction of the sky and will detect about N G = 5 x 10 7 
galaxies. We estimate the error in the power spectrum C\ measurement. The error in the 
IGM XRB-galaxy power spectrum is 



C G is the power spectrum of the galaxy surface density. C^(l) = AitffP ss /N G is the Poisson 



noise of the galaxy number count. C§(1) = 4nf^ ss (F N /F x ) 2 /N N + 4irf^t ss /N x , where 



the first term is the Poisson noise of the RASS background and the second term is the 
Poisson noise of the IGM XRB signal. JVjv and Nx are the total numbers of photons that 
the RASS received in the band 0.5-1 keV from background noise and the IGM, respectively. 
We adopt F N (0. 1-2.4 keV) ~5x 10~ 12 erg cm~ 2 s^deg' 2 (Fig. 3, Voges et al. (1999)). A/ 
is the bin width and jf^ ss reflects the cosmic variance. The strong non-Gaussianity of the 
XRB affects its error analysis. One could in practice estimate r)(l), the non-Gaussianity 
of the XRB-galaxy cross correlation power spectrum (for Gaussian case, rj — 0) from our 
models. For simplicity, we will adopt the result from our 512 3 moving mesh hydro (MMH) 
simulation (Zhang, Pen & Wang 2002). The parameters we adopted in our 512 3 simulation 
are fl = 0.37, Q\ = 0.63, VL B = 0.05, h = 0.7, a$ = 1.0, power spectrum index n — 1, box 
size L = lOO/i -1 Mpc and smallest grid spacing 40/i _1 kpc. Though the cosmology adopted 
in this simulation is not identical to the fiducial cosmology we adopt in this paper, it should 
yield a rough estimate of the XRB non-Gaussianity. During this adiabatic simulation we 
store 2D projection of the X-ray flux through the 3D box at every light crossing time through 
the box. The projections are made alternatively in the x, y, z directions to minimize the 
repetition of the same structures in the projection. Our 2D maps are stored on 2048 2 grids. 
After the simulation, we stack the XRB sectional maps separated by the width of simulation 
box, randomly choosing the center of each section and randomly rotating and flipping each 
section. The periodic boundary condition guarantees that there are no discontinuities in any 
of the maps. We then add these sections onto a map of constant angular size. Using different 




( V (l) + l)[Cf G ] 2 + (C? + C*) x (cp + Cg) 
(2/ + 1) AZ/gss 



(7) 



4 The factor of 4 is communicated to us by the referee Andrzej Soltan. 
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random seeds for the alignments and rotations, we make 40 maps of width 1.32 degrees by 
integrating from zero to z = 3 (Fig. 3). Using the same set of random seeds, we make 
40 maps of intergalactic gas surface density to measure the XRB-galaxy cross correlation 
assuming that galaxy number density traces the gas distribution. We then measure the 
fluctuations of the cross correlation power spectra of these maps to obtain the XRB-galaxy 
cross correlation power spectrum non-Gaussianity, namely, r)(l) (figure 4). The details of 
this simulation are described in Zhang, Pen & Wang (2002). The error is dominated by this 
non- Gaussian cosmic variance at large angular scales and by noise at small angular scales 
(fig. 5). C XG (l) can be measured to a better than 10% accuracy for 800 </ < 5000. If 
we cross correlate galaxies at redshift bin [zi — Az/2, Zi + Az/2] with the XRB and if Az is 
sufficiently small, C XG (l, Zl - Az/2,z, t + Az/2) -> P XG (J- y Zt ) x ^( Zi ) / A X (eqn. 2). This 
equation enables one to infer the redshift resolved Pxc{k,z). The error in the Pxc{k,z) 
estimation is given by eqn. (7) with all Ci being replaced by Cife — Az/2, Zi + Az/2) except 
for Cm, where should remain unchanged due to the absence of redshift information. 
We choose Az = 0.1. Pxc{k,z) can be extracted up to z ~ 0.6 and the error around the 
peak of A XG (k) = J^Pxaik) is < 20% (Fig. 5). Around this peak the dominant source 
of errors is noise and the effect of non-Gaussianity is negligible, the uncertainty of the rj(l) 
estimation does not affect the accuracy of Pxcik, z) around its peak and the subsequent 
<J 2 (z) extraction. The cross correlation coefficient r(k,z) = Pxc(k, z)/ \/ ' P x {k, z)Pa(k, z) 
has a weak dependence on z and enables one to infer the redshift dependence of Px{k,z) 
from the measurement of Pxg and Pq- Given this redshift dependence, one can invert the 
observable two-dimensional angular power spectrum C X {1) to a three-dimensional power 
spectrum Px(k,z). 

From these measurements, the feedback history can be extracted. The redshift averaged 
r c could be determined with a 10% accuracy (fig. 5). o"g as (z) oc Pxg( z ) could be determined 
with ~ 20% accuracy for z < 0.5 (fig. 5). Since ~ ^ *£* „ Msa, the feedback 

"gas ^gas T c %G 

amount E^q and the scale k gas or r c , at which feedback dominates over gravity, can in 
principle be extracted with a comparable accuracy. Its calibration would require simulations 
with feedback, which are currently being studied. 

Our estimation shows the sensitivity of the XRB statistics to the gas profile. In our 
estimates, the gas profile is taken as a free function and the gas fraction is taken to be the 
same for each halo. In practice, feedback changes both the gas profile and the gas fraction. 
This complication does not affect our estimation of the feasibility to extract the gas state 
from XRB observations, since it only depends on the sensitivity of the XRB statistics to 
the gas profile. But it does affect the interpretation of the data, for example, the relation 
between the XRB statistics and the amount of feedback. A further investigation of this issue 
requires a detailed study of the gas profile and its evolution when feedback presents. We 
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are currently carrying out a self-consistent calculation of the effect of the feedback on the 
XRB, the thermal Sunyaev Zel'dovich effect and the cluster L x -T relation and will check it 
in hydro simulations 5 . 

4. Conclusion 

We probed the feasibility to extract the IGM state and evolution from the combination 
of XRB surveys such as RASS and galaxy surveys such as SDSS. To do that, we first 
built two analytical models, the continuum field model and the halo model to calculate the 
statistics of the IGM XRB. The two models give consistent results on the IGM XRB flux 
and correlations. We found that the IGM may contribute a significant, if not dominant, 
fraction to the unresolved soft XRB flux, its auto correlation and cross correlation with 
galaxies. Since these statistics have different dependences on gas parameters such as the 
hot gas fraction and gas clumping factor, we suggest that by combining the XRB flux and 
correlation observations, hot gas fraction and hot gas clumping factor could be extracted 
simultaneously. At small scales, non-gravitational heating such as feedback from galaxies 
dominates over gravity. This changes the gas power spectrum and leaves signatures in the 
IGM XRB statistics and allows its extraction from XRB observations. From our models 
and the hydro simulation calibrated error estimation, we estimated that RASS+SDSS would 
constrain the gas clumping factor to a better than 20% accuracy up to z < 0.5. The amount 
of feedback and the scales where feedback dominates over gravity can be extracted with a 
comparable accuracy. 

Acknowledgments: We thank Andrzej Soltan and Wolfgang Voges for detailed discus- 
sions of the RASS. 
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